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Abstract 



We propose to modify the saturation model of Golec-Biernat and Wiisthoff by including 
DGLAP evolution. We find considerable improvement for the total deep inelastic cross 
section, in particular in the large Q 2 region. The successful description of DIS diffraction is 
preserved. 



The saturation model |T|, [2|, ||] has provided a successful description of HERA deep inelastic 
scattering data, in particular for the transition from the perturbative region to the nonper- 
turbative photoproduction region. This includes both the total j*p cross section and the DIS 
diffractive cross section. Whereas the formulae are particularly appealing through their simplic- 
ity, they also have an attractive theoretical background, namely the idea of saturation. Despite 
its success, the model suffers from shortcomings which should be cured. In particular, the model 
does not include logarithmic scaling violations, i.e. at larger values of Q 2 it does not exactly 
match with QCD evolution (DGLAP). This becomes clearly visible in the energy dependence of 
aj t in the region Q 2 > 20 GeV 2 where the model predictions are below the data. One expects 
that QCD evolution should enhance the cross section in this region. 

It is the purpose of this letter to propose a modification of the saturation model. We at- 
tempt to preserve the success of the model in the low-Q 2 and in the transition region, while 
incorporating DGLAP evolution in the large-Q 2 domain. Since the energy dependence in the 
large-Q 2 region is mainly due to the behaviour of the dipole cross section at small dipole sizes 
r, our changes will affect mostly the small-r region. At the same time, particular attention will 
be given to DIS diffraction for which the saturation model correctly describes the energy depen- 
dence. Since the inclusive diffractive cross section mostly depends upon the large-r behaviour of 
the dipole cross section, we attempt to leave the dipole cross section unchanged in this region. 
A recent attempt along the same lines indicates that, in fact, diffraction provides a highly 
nontrivial restriction on possible modifications of the saturation model. 



1 The Model 

Before we describe the modifications of the saturation model, we briefly review the main features 
of its original formulation. Within the dipole formulation of the 7*p scattering, 

a^l(x,Q 2 ) = J d 2 r J dz ip^ L (Q,r, z) a(x,r) *p T;L (Q,r,z), (1) 

where T, L denotes the virtual photon polarisation, the dipole cross section was proposed to 
have the form 

*(x,r) =a {l-exp(-^)}, (2) 

where Rq(x) is the saturation scale which decreases when x — > 0, 

1 / x x A 



EA(x) = ? — ■ (3) 

0V ' GeV 2 \x Q ) V ; 

In order to be able to study the formal photoproduction limit, the Bjorken variable x = xb was 
modified to be 

/ 4m 2 \ Q 2 + 4m 2 

where m q is an effective quark mass, and W 2 denotes the j*p center-of-mass energy squared. 
The parameters of the model, ao = 23 mb, A = 0.29 and xq = 3 • 10 -4 (for the assumed quark 
mass m q = 140 MeV) were found from a fit to small-x data Q. For alternative forms of the 
dipole cross section parameterisation see M . 
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As it is well known ||], in the small-r region the dipole cross section is related to the gluon 
density 

a(x,r) ~ — r a s xg(x,^ ). (5) 

where the scale fj, 2 for small r behaves as C/r 2 . Equation (5) is valid in the double logarithmic 
approximation in which the constant C is not determined. In the saturation model, eq. (g), we 
find for small r <C 2Rq{x) 

a(X ' r) " 4i^)' (6) 

i.e., the gluon density is modelled as: 

( 2\ 3 (T , . 

} = 4^o~ ' (7) 

For fixed a s this gluon density is clearly scale independent, which contradicts the QCD DGLAP 
evolution. Thus, in order to correctly take into account the scale dependence as given by the 
DGLAP evolution equations we have to modify the small-r behaviour of the dipole cross section 
by incorporating the properly evolved gluon density. At the same time, we wish to preserve the 
idea of saturation, which reflects unitarity, and to keep unaltered the large-r behaviour of the 
dipole cross section which determines the diffractive cross section. 

Therefore, we propose the following modification of the model (|2|) 

~r \ ft ( ir 2 r 2 a s (fi 2 )xg(x,fi 2 ) \\ 

a(x,r) = a < 1 - exp I — IS, (8) 

where the scale fi 2 is assumed to have the form 

/i 2 = ^ + /io- (9) 

The parameters C and /jlq will be determined from a fit to DIS data. In a first approximation, 
g(x,fi 2 ) is evolved with the leading order DGLAP evolution equation for the gluon density. In 
the spirit of the small- a; limit, we neglect quarks in the evolution equations. We assume the 
following gluon density at the initial scale Qq = 1 GeV 2 : 

xg(x,Q 2 ) = A g x-^{l-xf-\ (10) 

where A g and X g are parameters to be determined from a fit to data. The exponent 5.6, deter- 
mining the large x behaviour, is motivated by one of the versions of the MRST parameterisation 
|5J of the gluon density. 

For small r, the exponential in (§) can be expanded in powers of its argument, and the 
relation (||), with the running a s = a s ((j?), is found. In contrast to the original dipole cross 
section, the rise in 1/x now has become r-dependent. When inserting a into (|l]) and convoluting 
with the photon wave function, the integrand peaks near r ~ 2/Q for large Q 2 , and the argument 
of the gluon density turns into fi 2 ~ Q 2 . Consequently, with increasing Q 2 , DGLAP evolution 
will strengthen the rise in 1/x, whereas in the original saturation model the power of 1/x had 
been constant. For sufficiently large r, the scale fi 2 is frozen at the value /Uq. This prevents 
the effective scale of the gluon density from becoming unreasonably small. The saturation value 
of the dipole cross section is a(x,r) ~ uq, as in the original model (||). The transition from 
the small to the large r region depends on x, but in detail it will be different from the original 
model. This will be discussed in detail in the section presenting numerical results. 
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2 Momentum space formulation 



Although in this letter we will restrict ourselves to the total (and later on to the diffractive) 
cross section, it is instructive to rephrase these features in momentum space. In a future step, 
we intend to study the effects of saturation in more exclusive final states, and the translation of 
our dipole cross section into momentum space may serve as a first step into this direction. For 
this purpose let us start with the /c-r-fatorization formula ||] for the cross section, e.g for 
the transversely polarised photon M, 



a 



a ^r e e / r dz [z2 + (i - z)2] i d2k ii^ as/(x ' ^ 1 k k + 1 



" — 1 M ' " k 2 + Q 2 (k + l) 2 + Q 2 J ' 

(11) 

where f(x,l 2 ) is the gluon amplitude describing an interaction of the qq pair with the proton, 1 
is the transverse momentum of the gluon coupled to the quark pair and Q = z(l — z)Q 2 . Using 
the relation 

—^2 = iQ I ^ e ' k ' r £ Ki(Qr) , (12) 



k 2 + Q 2 J 2vr 

the following formula is found 

3a 



r 



a, 



7*p 
T 



2. 2 



^ e} £ dz[z 2 + (l- zf] Q 2 J d 2 r J d 2 r' ^ K t (Qr) K x (Qr') 



''" k ^< r -^D{r,r') (13) 



(27T) 



where 



D(r, = Y / -p- a./far, / 2 ) (1 - e* I r ) (1 - e ~< 1 " / ) . (14) 

If the argument of the strong coupling a s and the variable x in the gluon amplitude in D(r, r') 
do not depend on the quark transverse momenta k, the integration over k in eq. (|l^) gives the 
delta function 5 2 (r — r') which reflects the conservation of a dipole transverse size vector r during 
the collision. In this case the dipole formula (|T|) is obtained with the following identification of 
the dipole cross section 

a(x, r) = y / ^ a s f(x, I 2 ) (1 - e^) (1 - e" iI r ) . (15) 
Going beyond the leading log(l/x) approximation in which eq. ([ll]) was derived, e.g. by taking 



into account the exact gluon kinematics [1C] or considering a quark virtuality k 2 + Q 2 as an 



argument of the running coupling a s , we find that r is no longer conserved during the scattering 
process, and the simple relation (]i~5| ) ceases to exist. As a result, the /c^-factorization formula 
( |TT| ) can no longer be written in the form (|l]), and the simple dipole picture fails. We want to 
avoid this situation, thus we assume that the argument of a s is given by the gluon momentum 
l 2 and x = XBj- Since the integration in (|l~5|) includes also small momenta, the modelling of the 
infrared behaviour of a s cannot be avoided. However, we will hide this fact by analysing the 
combined quantity a s f(x,l 2 ). 



From the requirement that in the double logarithmic limit (DLL) formula (11) should be 



consistent with the DLL of the DGLAP formalism one can derive a relation between the gluon 
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amplitude f(x,l 2 ) and the conventional gluon distribution xg(x,Q 2 ). Starting from eq. (|Tl|), 
using the relation Ft = Q 2 / (4:ir 2 a em ) a T p and imposing the strong ordering condition: l 2 <C 
k 2 <C Q 2 , one arrives at 

dF T {x,Q 2 ) _ 1 2 rQ 2 d 2 \ 



^5>?/ ^sfixA. (16) 



<91ogQ 2 ,,. , 

By comparison with an analogous formula obtained in the DLL of the DGLAP evolution equa- 
tions, one finds the following relation at large Q 2 

a s (Q 2 )xg(x,Q 2 ) = ^ fia s f(x,l 2 ). (17) 

In the model @ we go beyond the fc^-factorization formula ( |Tl|) where f(x,l 2 ) represent a 
two gluon amplitude. In the region of small I 2 , the relation ( |i~7|) between f(x, I 2 ) and the gluon 
density no longer holds, and f(x,l 2 ) is defined through relation ( [l5|) where for the l.h.s we use 
our model (§). In general, provided the dipole cross section has a finite limit: lim^oo a(x,r) = 
&oo( x )i t ne equation (FL5|) can be inverted with the help of the following relation: 



a s f(x,l 2 ) 3fd 2 



f d r 



l A 4vr J (2tt) s 

drrJ (lr) {&oo(x) - a(x,r)} . (18) 



_3_ 

8^ 2 Jo 

In the original dipole model we find Q: 

a s f(x, I 2 ) = R 2 (x) I 4 eM-Rfa) I 2 } ■ (19) 

For the modified dipole cross section, this inversion has to be done numerically. DGLAP evo- 
lution will affect mainly the large-/ behaviour while at small I our modification should be less 
severe. The most interesting question to be addressed below concerns the transition region: to 
what extent does our modification affect the region of moderate momenta, i.e. could one 'see' 
saturation in diffractive final states ? 



3 Numerical results 

Let us now turn to numerical results. We performed global fits to the DIS data with x < 0.01 
in the range of Q 2 between 0.1 and 500 GeV 2 . For HI and ZEUS HERA experiments the new 



1996-1997 data sets were used [11, 12, [DJ. In addition to the HERA data also the data of the 
E665 experiment H] were used. The statistical and systematic errors were added in quadrature. 
The number of degrees of freedom, Njf, was around 330. 

The new data sets are considerably more precise (with much smaller statistical and systematic 
errors) then the ones used in the original analysis As a preparatory step, we applied the 
original model @, using the parameter values of the original fit, to the new data and obtained 
rather high value of X 2 /N^f ~ 3 (for the old data, the corresponding value was X 2 /N^f = 1.18). 
Next, we let the new data to determine their own values of the parameters of the original model, 
(Tq, A and xq in eq. (||). This led to an improvement of the fit, X 2 l^df ~ 2.2. Nevertheless, this 
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m q (MeV) 


o"0 (mb) 


A 9 


\ 


C 




X 2 /N df 


Fit 1 


140 


23.0 


1.20 


0.28 


0.26 


0.52 


1.17 


Fit 2 





23.8 


13.71 


- 0.41 


11.10 


1.00 


0.97 



Table 1: The parameters of the fits to the ZEUS, HI and E665 data with x < 0.01 (333 points). 
The HI data was rescaled by a factor of 1.05. The numbers in bold are fixed during the fits. 

relatively poor agreement indicates that the original model is doing not so well with the new 
data, especially for large values of Q 2 . As a first step for the improvement, we modify the dipole 
cross section at small values of r by including QCD DGLAP evolution, as given in eq. (Q). 

In the modified saturation model, there are five parameters to be determined: do, C, (i 2 ,, A g 
and X g from eqs. (|8|, ||, [R]). We use the leading order DGLAP evolution equation for the gluon 
density, and we put Aqcd = 200 MeV in a s and set the number of active flavours Nf = 3. Thus, 
although the evolution equation for the gluon is decoupled from the quarks, their presence is 
encoded in the assumed value of Nf. 

We performed first the fit leaving all five parameters free and assumed the value of the light 
quark mass m q = 140 MeV, as in the original formulation GJ. A good quality fit was obtained 
with x 2 /Ndf = 1.05. The found value of the dipole cross section <to = 27.4 mb, however, was 
higher than the saturation model value, do = 23 mb. Also, the corresponding value of the 
photoproduction cross section (cr' rP = 204 mb) was significantly higher than the measured value 
174 ± l(st) ± 13(sys) mb at W = 209 GeV p|. Thus we decided to decrease a and fixed it 
to the saturation model value 23 mb. This value is also advantageous for the description of the 
inclusive DIS diffractive cross section which is more sensitive to large dipole sizes, i.e., to the 
saturation region, than the total 7*p cross section PJ. The resulting parameters in such fit are 
presented in Table 1 (Fit 1). The description is slightly worse than that described above, but 
both the photoproduction cross section (<t 7P = 189 mb) and the diffractive cross section are 
properly described. This is because we modified only the small dipole size part of the dipole 
cross section (g), without affecting the saturation part. The found gluon density gives 39% of 
the total proton momentum carried by gluons resolved at the initial scale Qq = 1 GeV 2 . 

The results of Fit 1 are compared to the data on F2 in Fig. |l| for Q 2 < 1 GeV 2 and in 
Fig. ^ for large Q 2 points. In all presented plots, the solid lines refer to the results obtained 
with the DGLAP improved model (|8|) and the dashed lines correspond to the saturation model 
(H|) with the original parameters from We see that the DGLAP evolution significantly 
improves agreement with the data at large Q 2 while at small Q 2 the results are practically the 
same. This effect is summarised in Fig. [3| where the effective slopes \(Q 2 ), obtained from the 
parameterisation of F2 at small x: F2 ~ x~ x ^ \ are plotted. Thus, the DGLAP modification of 
the dipole cross section for small r is crucial for much better agreement with the data. The same 
effective slopes characterise the energy dependence of the -y*p cross section: cj 7 * p ~ (W 2 ) x ^ >. 
The change from a soft dependence at small Q 2 to a hard one for large Q 2 is shown in Fig. [|[ In 
Fig. [| we show another aspect of the transition of F2 to low Q 2 values, namely the emergence of 
the behaviour: F2 ~ Q 2 approached in the limit Q 2 — > and y = W 2 / s fixed. The class of the 
saturation models described here nicely reproduce this behaviour, see recent Ref. |l6|] for more 
details on this transition. 

In a second step of our investigations we relax our requirement of staying in the low-Q 2 
region as close as possible to the original model. In particular, we allow the quarks in the qq 
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dipole to become massless. Thus we set m q = in the wave function ^>t,l and in the kinematic 
relation (|j). In the original model, the quark mass was introduced as an effective parameter for 
modelling the large-r behaviour of the photon wave function. The non-zero quark mass allows 
to study the photoproduction limit of our model after the modification (Q) of the Bjorken-x in 
the dipole cross section, therefore, setting this parameter to zero eliminates this possibility. It 
allows, however, for a better description of the current data. We also fix the minimal value \£q 
of the scale /i 2 in eq. ((^) to 1 GeV 2 in the fits in order to avoid negative gluon density below 
the input scale Q\ = 1 GeV 2 for the gluon evolution^. Since the parameters gq and X g are 
strongly correlated, we have performed a systematic search of the best % 2 on the grid of fixed 
(do, X g ). In each case, the remaining parameters, A g and c, were fitted. In this way we found 
two local minima for % 2 , shown in Fig. ||, for positive value of X g in eq. (fL0|) leading to strongly 
rising gluon density, and for negative value of X g , corresponding to valence- like initial gluon. 
The latter scenario gives considerable better description (with x 2 /N^f = 0.97) than the first one 
(with x 2 /Ndf = 1.13). In the final analysis, after a quantitative estimation of the position of the 
best fit in the parameter space using the grid method, we allow <jq and X g to be fitted together 
with A g and c. The values of these parameters for the best fit are given in Table 1 (Fit 2). The 
corresponding value of the gluon momentum at the input scale equals 84%. The effective slope 
X(Q 2 ) from the Fit 2 parameterisation is shown as the dotted line in Fig. ||[ As expected, slight 
differences between the two fit scenarios only appear for small values of Q 2 , below 1 GeV 2 . 

It is interesting to compare the results of the Fit 1 and Fit 2 since they lead to a different 
picture of the dynamics of the "f*p interaction. In Fit 1, the initial gluon density, xg(x,QQ), 
quickly rises with 1/x, whereas in Fit 2 it even decreases with rising 1/x. Therefore, in Fit 1 
the rise of the cross section with the energy is mainly due to the intrinsic properties of the initial 
gluon density, with only slight corrections being due to the evolution effects at high Q 2 values, 
and considerable damping effects resulting from saturation at low Q 2 . In Fit 2, the evolution 
effects are very strong (note the value of the parameter C, which is much higher then in Fit 1). 
The small- x rise of the cross section is due solely to the DGLAP evolution effects with some 
corrections coming from saturation. 

Further insight into the physical picture lying behind the fits can be gained by a closer look 
at r-dependence of the dipole cross section a(x,r), and a momentum dependence of the related 
gluon amplitude, f(x, I 2 ). In the first row of Fig. ^ we show the dipole cross section for the two 
fits. As we have already discussed, for Fit 1 the DGLAP modification (solid lines) affects mostly 
the region of moderately small r (< 1 GeV -1 ). In Fit 2 both the small and large r regions are 
affected. In particular, the structure close to the saturation region is different from that in the 
saturation model (|2|) (dashed lines). The differences between the models are particular visible 
if we turn to momentum space and compute the gluon amplitude a s f(x, I 2 ) for different values 
of x, using relation (|l8|). The results are shown in the second row of Fig. [?], where the full lines 
denote the gluon amplitude from the DGLAP improved model and the dashed lines correspond 
to the saturation model (|2|). The small r region of the dipole cross section corresponds to the 
large I 2 region in the gluon amplitude. The dipole cross section from Fit 1 is translated into 
a gluon amplitude with a double-bump structure. Notice that the second bump results from 
the DGLAP modification of the small-r part of the dipole cross section. In Fit 2, however, the 
second bump disappears and the gluon amplitude is similar but significantly broader than the 
one corresponding to the saturation model. Although in Fig. [7| the various gluon amplitudes are 
clearly distinct, after convolution with the impact factors and turning to the 7*p cross sections, 

The valence-like gluon density preferred in the massless fit, as described in the text, becomes negative below 
the input scale due to backward evolution. In this case the dipole cross section (feh does not saturate at large r. 
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these differences are becoming much less visible. It is natural to expect, that some of these 
differences should become visible in more exclusive final states. As a first example, one might 
think of DIS diffraction. The inclusive diffractive process, however, is sensitive to the region of 
small I 2 or large r and, therefore, has only limited value in distinguishing between the two fit 
solutions. On the other hand, other physics processes like jet, charm and bottom production 
should be more sensitive to the behaviour of the unintegrated gluon density at large gluon 
momenta I 2 . 

Looking at Fig. fj] notice that starting from certain values of I 2 , the gluon amplitudes become 
negative. In order to understand this let us differentiate the relation (|l7]) with respect to the 
large scale fj, 2 , 

a ^" 2 > = ah W" 2 >*" 2 >} " a '<" 2 > ~ 00)} ■ < 2 »> 

The quantity in the curly brackets in the last equality can become negative, which is shown 
in the bottom row of Fig. [?] by plotting the r.h.s of the above equation (dashed lines). In the 
shown range of I 2 , eq. (|20| ) is especially well satisfied for the parameters from Fit 2. For Fit 1 
the equality is reached for much larger (not shown) values of I 2 = fj, 2 . 

In ref. Q the critical line in the (x, Q 2 ) plane was defined which marks the transition to the 
saturation region where a new behaviour of the structure function, F2 ~ Q 2 , emerges. Near this 
line, the characteristic size of the qq dipole, f ~ 2/Q, equals the saturation radius Rq(x), see 
section |l[ In this case the argument of the exponent in eq. (||) equals one, 

Q 2 R 2 (x) = 1 (21) 

and o"(x,f) ~ Co- We adopt the same criterion for the critical line in the DGLAP improved 
saturation model (@). Thus, we have the following condition 

47T 2 

a s {ii 2 )xg{x^ 2 ) = 1, (22) 



3 do Q 

where fi 2 = C Q 2 /4 + //q. Eq. (22) is an implicit equation for the critical line Q 2 = Q 2 (x 



shown in Fig. |8] for the two fits. As expected, for FIT 1 the found critical line is not different 
from that defined in the original saturation model (dashed line) and the transition region stays 
around 1 GeV 2 in the HERA kinematics (lower band). For FIT 2 the critical line is situated at 
lower values of Q 2 (around 0.5 GeV 2 ). It is interesting to note that both fits predict that in the 
THERA kinematic range (upper band) the saturation region lies at Q 2 ~ 2 GeV 2 , which puts 
the perturbative QCD description of saturation effects on more solid ground. 



4 Diffraction 

One of the main advantages of dipole models is their straightforward description of diffractive 
processes. The generalised optical theorem applied in the framework of the dipole picture allows 
to express the cross section for diffractive qq production in which proton remains intact as 



dt 1 t=o 



— J d 2 r J dzr T ,dQ,r,z)a 2 {x,r)^ L {Q : r : z), (23) 
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where t = A 2 , and A is the four-momentum transfered into the diffractive system from the pro- 
ton. In addition to the contributions of the qq states it is important to include the contributions 
of the qqg final states [17]. 

In the phenomenological analysis 0] , the qqg diffractive amplitude was computed in the two- 
gluon exchange approximation with an additional assumption of strong ordering of transverse 
momenta of the qq pair and the gluon. This allows to treat the qqg system as a color octet 
dipole (88) in the transverse coordinate representation. Compared to the triplet dipole, the 
coupling of two t-channel gluons in the singlet state to the octet dipole carries the relative 
weight Ca/Cf = 2N 2 /(N 2 — 1). Thus, in order to take into account the repeated exchange of a 
two gluon system, the equation (||) for the triplet dipole cross section is modified for the octet 
dipole as 

a gg (x,r) = a jl - exp — — j j . (24) 

The above modification is done in the spirit of multiple pomeron exchange, i.e. the term pro- 
portional to (r 2 ) n , resulting from the expansion of the exponent in (p4|), would correspond to 
n exchanged pomerons with an appropriate colour factor (Ca/Cf) 71 ■ In addition, compared to 
the diffractive qq production, the cross section formula for a diffractive qqg system contains an 
overal factor (N 2 -l)/N c . 

One of the most important results of the original saturation model was that, at fixed Q 2 , the 
ratio of the inclusive diffractive cross section and the total j*p cross section is nearly constant in 
agreement with data. This prediction is not changed in the DGLAP improved saturation model 
since we modified only the short distance part of the dipole cross section. Even in the case of 
Fit 2, the constant ratio is preserved, as shown in Fig. [9|, in contrast to the attempts in Q. The 
theoretical curves in these figures are computed using the Fit 2 results, and the experimental 



data are taken from ZEUS [18]. The results for the Fit 1 computation differ only slightly from 
the Fit 2 one. 

The diffractive data shown in Fig. |9| were obtained without the experimental identification 
of the forward going proton. Therefore, as described in [18], this data have a substantial contri- 



bution of the proton dissociation process which was estimated as 31 db 15%. To take into account 
this contribution, the prediction of our model shown in Fig. ^ were multiplied by a normalisation 
factor of 1/(1 — 0.31) = 1.45. The agreement of the hight of the predicted cross section with 
the data is satifactory only within the relativly large error of the estimated proton dissociation 
factor. 

We also made a comparison of our predictions with the recent preliminary diffractive data in 
which the forward going proton was identified in the ZEUS Leading Proton Spectrometer (LPS) 
[19]. Good agreement was found which gives further support to the dynamical picture of the 



7*p interactions developed in this and previous papers on the saturation model. 



5 Conclusions 

In this letter we have proposed a modification of the saturation model which takes into account 
the QCD DGLAP evolution of the gluon distribution. Fitting the parameters of our model we 
found a solution that describes the new HERA data on F% significantly better than the original 
saturation model, especially in the region of larger Q 2 . The agreement with the DIS diffractive 
HERA data is also kept. 
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Somewhat surprisingly, we found another set of parameters which lead to even better data 
description departing from the original saturation model. For this description, we set the effective 
quark mass of the original model equal to zero, and in our comparison with the HERA data 
we have disregarded photoproduction data points. We found indications that this solution 
represents a slightly different physical picture: the initial gluon density no longer rises at small 
x, and QCD-evolution plays a much more significant role than in our first solution. The fact 
that the effective quark mass of the original model has its strongest influence in the limit Q 2 — ► 
suggests that the large-r behavior of the photon wave function requires further considerations. 

We have found it useful to discuss the various versions of the saturation model not only in 
the r-space but also in momentum space since the latter provides more direct connection with 
exclusive final states. As a future step, it will be instructive to trace saturation effects in less 
inclusive cross sections. 

We consider the modification of the saturation model presented in this paper as a first step of 
a more systematic program. The success of the original model indicates that this simple ansatz 
contains elements of the correct dynamics. Next, we have to analyse this model within QCD 
and to find the necessary corrections. With precise HERA data on various reactions becoming 
available available, all modifications have to be testet by careful comparisons. 
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ZEUS BPT97 




Figure 1: F% as a function of x for fixed low Q 2 values. The comparison with the low Q 2 data 
from ZEUS. The solid lines: the model with the DGLAP evolution ^) (FIT 1) and the dotted 
lines: the saturation model 
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Figure 2: HI and ZEUS data on F2 as a function of x for fixed values of Q 2 > 1 GeV 2 and the 
saturatiom model curves. The solid lines: the model with the DGLAP evolution (^) (FIT 1) and 
the dotted lines: the saturation model 
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Figure 3: The effective slope X(Q 2 ) from the parameterization F2 ~ x~ A (<2 2 ) as a function of Q 2 . 
The model with the DGLAP evolution (^): the solid line (FIT 1) and the dotted line (FIT 2). 
The saturation model (^): the dashed line. The open circles: ZEUS analysis and the full circles: 
HI data JMj. 
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Figure 4: The j*p cross section as a function of energy W 2 at various Q 2 . The solid lines: the 
model with the DGLAP evolution (FIT 1 ) and the dotted line: the saturation model (||), 
shown for x < 0.01. 



14 



Transition to low Q 2 



CM 

o 

W 

CM 

o 

II 

^10 3 



10 



10 



10 



y=0.8 
(x 4096) 

y=0.7 
(x 2048) 

y=0.6 
(X1024) 

y=0.5 
(x 512) 

y=0.4 

(x 256) 

y=0.33 
(x 128) 

y=0.26 
(x 64) 

y=0.20 

(x 32) 

y=0.12 
(X16) 

y=0.05 

(x 8) 

y =0.025 
(x 4) 

y=0.015 

(x 2) 

y =0.007 
(x1) 




10 



10 



10 



10 

Q 2 (GeV 2 ) 



Figure 5: F2(x,Q 2 ) as a function of Q 2 for fixed y = Q 2 /{sx). The solid lines: the model with 
DGLAP evolution (FIT 1) and the dashed lines: the saturation model (j^J. The curves are 
plotted for x < 0.01. Full circles: ZEUS data and open circles: HI data. 
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Figure 6: The lines of constant values ofx 2 in the space of (<ro, — X g ) in the massless case, m q = 0. 
The two local minima are indicated by the black dots. The one for —\ g — 0.4 corresponds to 
FIT 2. 
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Figure 7: The dipole cross section a(x,r) (upper row) for x = 10~ 2 ... 10~ 6 (from right to left) 
and the gluon amplitude a s f(x,l 2 ) at x = 10~ 2 ... 10~ 4 (from bottom to top) for the two fits. 
The solid lines corespond to the DGLAP improved model while the dashed lines describe the 
saturation model (Q). The dotted lines in the bottow row show the r.h.s of eq. (|2C|). 
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Figure 8: The position of the critical line in the (x, Q 2 ) plane in the DGLAP improved model 
(solid lines) and the original saturation model (dshed line). The bands indicate acceptance 
regions for the colliders HERA (lower) and future THERA (upper). 
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Figure 9: The ratio of (Jdiff/^tot versus the 7*p energy W. The data is from ZEUS and the solid 
lines correspond to the results of the DGLAP improved model with massless quarks (FIT 2). 
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